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Output Feedback Control of the Linear Korteweg-de Vries Equation* 


Swann Marx 1 and Eduardo Cerpa 2 


Abstract — This paper presents the design of an output 
feedback control for a linear Korteweg-de Vries equation. This 
design is based on the backstepping method which uses a 
Volterra transformation. An appropriate observer is introduced 
and the exponential stability of the closed-loop system is proven. 

I. Introduction 

The Korteweg-de Vries equation, introduced in 1895 by 
the Dutch mathematicians Diederik J. Korteweg and his stu¬ 
dent Gustav de Vries, describes approximatively the behavior 
of long waves in a water channel of relatively shallow depth. 
This nonlinear partial differential equation, described by 

Ut(t, X ) + U x (t, X ) + u xxx (t, x) + u{t, x)u x {t , x) = 0, (1) 

has been deeply studied in the controllability sense (see [9], 
[4], [1] and the references therein). By considering diffe¬ 
rent boundary conditions on an interval [0,L] and different 
boundary actuators, we get control results of different nature. 
Roughly speaking, the system is exactly controllable when 
the control acts from the right endpoint x = L, and null- 
controllable when the control acts from the left endpoint x = 
0. 

In this article we focus, as a first step opening further re¬ 
search, on the stabilizability problem for the linear Korteweg- 
de Vries equation with a control acting on the left Dirichlet 
boundary condition. The studied system can be written as 
follows: 

tit + ti x T u xxx — 0, 

u(t, 0) = k( t), u(t , L ) = 0, u x (t, L) = 0, (2) 

w(0,x) = u 0 {x), 

where k denotes the control input and u o is the initial 
condition. 

Some full state feedback controls have already been de¬ 
signed in the literature. Let us mention [3] where a Gramian- 
based method is applied in the case where the control acts on 
the right endpoint and [11] where the backstepping method 
is applied with different boundary conditions. With the same 
approach than the last paper, we have [2] where system ([2} 
is considered. 

However, in most cases, we have no access to the full 
state of the system, and it is more realistic to design an 

*This work has been partially supported by FONDECYT grant 1140741, 
MathAmSud COSIP, and CONICYT grant ACT-1106 

1 Swann Marx is a student from the Electrical Engineering Department, 
Ecole Normale Superieure de Cachan, 61 Avenue du President Wilson, 
94230 Cachan. France swann . marxSens-cachan . fr. 

2 Eduardo Cerpa is with Departamento de Matematica, Universidad 
Tecnica Federico Santa Marla, Avda. Espana 1680, Valparaiso, Chile 
eduardo.cerpa@usm.cl. 


output feedback control, i.e., a feedback law depending only 
on some partial measurements of the state. 

For linear and autonomous finite-dimension systems, a 
stabilizability and an observability assumptions are sufficient 
to ensure that the separation principle holds. In other words, 
if there exists a controller, which asymptotically stabilizes 
the origin of the system and an observer which converges 
asymptotically to the state system, the output feedback built 
from this observer and this state feedback asymptotically 
stabilizes the origin of the system. In a PDE framework this 
principle is no longer true and the stability of the closed-loop 
system is not guaranteed. 

The basic question to state the problem is which kind 
of measurements are we going to consider. The case of a 
boundary measurement is the most challenging one and the 
natural choice for the KdV equation ({2} should be y(t) = 
u x {t,0). Unfortunately, the system is not observable with 
this choice. In fact, if L = 2n, then the stationary solution 
u(t,x) = 1 — cos(x), satisfies 0 with n(t) = 0 and in 
addition u x (t, 2it) = 0 for any time t. The length 27r is not 
the only one for which observability does not hold (see e.g. 
[8], [1]). 

In this paper we consider the output given by 

y{t) =u xx (t,L). (3) 

By using this measurement, we build an observer and apply 
the backstepping method to design an output feedback con¬ 
trol which exponentially stabilizes the closed-loop system. 
The trace u xx (t,L) can be considered itself as a boundary 
condition or even a control ([5], [4]) and therefore the choice 
of this output is not an artificial one. 

This paper is organized as follows. In Section [HI we 
state our main result. Section [HI] is devoted to recall the 
state feedback control designed in [2], In llVI we state some 
regularity results needed to consider the measurement (0 
as a continuous function. The observer is built in Section 
IV1 In Section [VI] the stability of the closed loop controller- 
observer system is proven. Finally, Section IVIII states some 
conclusions. 

II. Main Result 

Based on [7] and [10], we construct the following ob¬ 
server: 

tit ~b tix T Uxxx T Pi (*r) [y(f) ti xx (t 1 Z/)] — 0, 
u{t, 0) = n(t), u(t, L) = u x (t , L) = 0, (4) 

u(0, x) = 0. 

Theorem 1. For any A > 0, there exist a feedback law 
n(t) := n(u(t,x)), a function pi = p\(x), and a constant 


C > 0 such that the coupled system ©-©-© is globally 
exponentially stable with a decay rate equals to A, i.e., for 
any uq £ H 3 (0 : L) we have 

ll u (*> •)llff 3 (0,i) + l|w(^ OIU^O.L) <Ce At ||uo||_H'3( 0 i ) 

(5) 


III. Full state feedback design 


In [2] the following Volterra transformation is introduced 

[ L 

w(x) = n(rt(a;)) := u(x) — / k(x,y)u(y)dy. (6) 

J X 

The function k is chosen such that the trajectory u = u(t, x), 
solution of ©, with control 

n{t)= [ k(Q,y)u(t,y)dy, (7) 

Jo 

is mapped into the trajectory w = w(t,x), solution of the 
linear system 


Wt + Wx + Wxxx + \w = 0 , 

w(t, 0) = w(t, L) = w x (t, L) = 0, 


( 8 ) 


which is exponentially stable if A > 0. In fact, from a 
Lyapunov approach, we get 


d_ 

dt 



| w(t, x)\ 2 dx 


= —\w x {t, 0)| 2 — 2A 



\w(t, x)\ 2 dx 


< — 2A f \w(t,x)\ 2 dx, 

Jo 


(9) 


which gives an exponential decay rate equals to A for the 
L 2 — norm of the state w. 

The kernel function k = k(x, y) is characterized by: 

k X xx 4- kyyy 4- k x T ky = A/c, in T , 
k(x, L) = 0, in [0, L\, 

* k(x,x) = 0, in [0,1/], 

kx{x,x) = ^(L - x), in [0,L], 

where T := {(x,y)/x £ [0 ,L\,y £ [x,L]}. The solution of 
(US exist. This is proved by using the method of successive 
approximations. Unlikely the case of heat or wave equations, 
we do not have an explicit solution. 

In [2] it is proved that the transformation © linking 
(0 and 0 is invertible, continuous and its inverse is also 
continuous. Therefore, the exponential decay for w, solution 
of ©, implies the exponential decay for the solution u 
controlled by 0. Thus, with this method, the following 
theorem is proven. 

Theorem 2. (State feedback stabilization for KdV ([2]) 

For any A > 0, there exist a feedback control law n = 
ti(u(t, .)) and C > 0 such that 

I \ u ( t , OHl^o.l) < Ce~ xt ||uo||l 2 (o,l) (H) 

for any solution of 


IV. Regularity Result 
A s we said in the introduction, we consider 

y(t) = u xx (t, L) (12) 

as a partial measurement of the solution. However, since we 
have the trace of the second derivative with respect to x of 
u, we need a regularity stronger than in [2], Indeed, we ask 
the output y(t) to be a continuous function. Thus we have 
the following lemma. 

Lemma 1. Let us consider system 

{ lit T ti x u xxx = 0, 

u(t, 0) = K (f), u(t, L ) = 0, u x (t, L ) = 0, (13) 

m(0,x) = u 0 (a;), 

where Uq £ 7T 3 (0, L) and n(t) £ fT 1 (0,T). Then u £ 
C([0, T], 7T 3 (0, L)) n L 2 (0,T;ff 4 (0,L)) and u xx {-,L ) £ 
C([0,T]). 

Proof This proof is based on [6]. Let us consider the 
following coordinates transformation 

v = u t . (14) 

The dynamics of v can be written as follows: 

{ rtt T v xxx -\- v x = 0, 

v{t, 0) = k(t) £ L 2 (0,T), v(t,L) = v x (t,L) = 0, 
v(0,x) = (-u"- u' 0 )£L 2 (0,L). 

(15) 

By already known well-posedness results for KdV ([6]), 
we get 

v £ C([0, T\,L 2 { 0, L)) n L 2 (0, T; H\0, L)) 

^u t £ C([0, T], L 2 (0, L)) n L 2 (0, T; H\0, L)) 

^u £ C([0, T\,H 3 ( 0, L)) n L 2 (0, T; H\ 0, L)) 

i o o (lo) 

^u xx £ C([0, T], H\0, L)) n L 2 (0, T; H 2 ( 0, L)) 

=^u xx £ C([0, T] x [0, L\) 

=► U XX (;L)£C([0,T ]) 

Thus it concludes the proof of Lemma Q] □ 

V. Observer design 

Based on [7], and more precisely on [10], we can write, 
for system 

Ut T u x + u xxx = 0, 

< u(t, 0) = n(t), u(t, L) = u x (t, L) = 0, (17) 

. y{t) = u xx (t,L), 

the corresponding observer 

j Ut 4“ u x 4~ u xxx 4~ Pi (a:) [rta;a: (U L') u xx (t , L)] — 0, 

1 u(t, 0) = «(f), u(t, L) = u x (t, L) = 0. 

(18) 

The construction of the observer is based on the finite- 
dimensional design for x = Ax + Bu, y = Cx, which 
proposes the observer x = Ax 4- Bu 4- L(y — Cx). If 
we consider the error e = x — x, then e = {A — LC)e, 




By adding (l22l >. (f23l) and (l25l >. we get 


and we have to look for a matrix L insuring a good 
performance. Because of the infinite-dimensional framework 

we are working in, a matrix L is not enough. Thus we need u t +u x + u xxx — Pi{x)u xx (L) — 
a function pi(x). 

In our case, we consider the error u := u — u, which 
satisfies 


w t (t,x) + w x (t,x) + w xxx (t,x) + \w(t,x) 


(19) 


j U t +u x + Uxxx - Pl{x)u XX (t, L) = 0, 

| u(t, 0) = u(t, L) = u x (t, L) = 0. 

Given a positive parameter A, we look for a transformation 
n o defined by 

u{x) = n o (w(a;)) = ui(x) - / p{x,y)w(y)dy (20) 

J X 

such that the trajectory u, solution of ( fl9l > is mapped into 
the trajectory w = w(t,x), solution of the linear system 


- / (~Xp{x,y) +p v {x,y) +p vvy (x,y) 
J X 

+ Pxxx{x,y) +Px(x,y))w(y,t)dy 


w t +w x + w xxx + Xw = 0, 

w(t, 0) = 0, w(t, L) = 0, w x (t, L) = 0, 


( 21 ) 


which is exponentially stable with a decay rate depending 
on the value of A as shown in ©. 

Now, the key step is to find the kernel p = p(x, y) such 
that u(t, x) = U 0 (w(t, x)) satisfies (fl9l l. By focusing on (f20l) 
and using the Leibniz rules, we get: 

• Differentiation along (fTil) 

u t = w t {x)~ / p(x,y)[-w y (y,t)-wyyy{y,t)-\w{y,t)\dy 
J X 

= w t (x)~ / (—Xp(x, y) + p y (x, y) + p yyy (x, y))w(y, t)dy 
J X 

+ p(x, L)w(L , t) - p(x, x)w(x, t) + p(x, L)w xx (L, t) 

-p(x,x)w xx (x,t ) +p y (x,x)w x (x) -p y (x,L)w x (L,t) 

+ Pyy(X,L)w(L,t ) ~ Pyy(X,X)w(X,t) 

( 22 ) 


+ w x (t, x) \ 2—p(x, x) + p x {x, x) + Py(x, x) 

+ w(t, x) ^p xx (x, x) + -^p{X, x) + -^p x (x, x) 

-p yy {x,x) - A) +p(x,L)w(L,t) + ( p{x,L) - Pi{x))w xx (L,t) 
-p y (x,L)w x (L,t). 

(26) 

From this equation, we get four conditions: 

1. Equation for (x,y) G T: 

Pyyy{x,y) +Pxxx(x,y) +p y {x,y) +p x (x,y) = \p{x,y). 

(27) 


2. First boundary condition on (x, x) for x £ M: 


2 —p(x,x) +p x (x,x) +p y (x,x) = 0. 


(28) 


3. Second boundary condition on (x,x) for x G [0, L\. 


_p( x , x) + _p, < x , x ) 

+ Pxx{x, X) ~ Pyy(x, X) - A = 0. 

4. Appropriate choice of pi. 

p{x, L) =pi(x). 


(29) 


(30) 


• Three differentations with respect to the variable x 

u x (x,t) = w x (x,t)+p(x,x)u>(x,t) - / p x (x,y)w(t,y)dy 

J X 

(23) 


U xx (t, x) = w xx (t, x) + —p{x, x)w(t, x) + p(x, x)w x {t, x) 
dx 


Recall that T := {(x,y)/x G [0, L], y G [x,L]}. 

Moreover, note also that, by setting x = 0 in (l20l) . we get: 


+ p x (x,x)w{t,x) - / p xx (x,y)w(t,y)dy 
J X 


(24) 


U xxx (t,x) = w xxx (t,x) + -^p(x,x)w(t,x) 

+ 2 -^-p(x,x)w x (t,x) 
ax 

+ p(x,x)w xx (t,x) + —p x (x,x)w(t,x) 
+ p x {x,x)u> x (t,x) +p xx (x,x)w{t,x) 


p(0, y) = 0, My G [0, L]. 

In addition, we have 

w(t, L ) = u(t, L) = w x (t , L) = 0. 

Finally, the kernel p satisfies the following PDE: 

' Pxxx(x,y) +Pyyy{x,y ) 

+ Py{x, y) + Px(x, y) = A p(x, y), (x, y) G T , 

p(x, x) = 0, x G [0, L\, 

p x (x,x) = ^x, xg[0,L\, 

,F(0,y) = 0, yG[0,L]. 

Let us make the following change of variable: 


(31) 


(32) 


- / p xxx (x,y)w(t,x)dy 

J X 


x = L-y, y = L-x, 


(33) 




and define F(x,y) \=p(x,y). Hence: 

F XXX (•£; y) T F yyy (X, ?/) 

+ F v (x, y) + Fx(x, y) = -A F(x, y) (x, y) € T 
< F(x,x) = 0 i 6 [ 0 ,L] 

Fx(x,x) = -x) xe[0,L\ 

,F(x,L) = 0 ye[0,L] 


Remark 1. We can prove that this Lyapunov function is 
equivalent to the one obtained by replacing Vs (t) by 


v 3 (t) 


B 

"2 





(t, x )\ 2 dx 


(41) 


By this, we mean that the exponential decay of one of them 
implies the exponential decay of the other one. In fact, we 
can prove that there exist positive constants d\, c ?2 such that 


This PDE has already been studied in [2], where no 
explicit solution has been found, but where the existence of 
a solution has been proved. Hence, we can conclude that the 
kernel p := p(x, y) exists. Note that the function n o defined 
by (l20l) is linear (by definition) and continuous (because of 
the existence of p). 

VI. Stability analysis of the closed loop system 

Instead of dealing directly with the controlled state u and 
the observer state u, we consider the evolution of the couple 
(u, u) where u stands for the error u = u — u, as introduced 
in Section [V] 

By using the output feedback control 

«(*) = [ ,y)u{t,y)dy, (35) 

Jo 

the transformation n defined in ([6]i and its inverse, and the 
transformation H Q defined in (l20l > and its inverse, we can 
see that (u,u) are mapped into (w,w) = (n~ 1 ({i), n(u)) 
solutions of the target system 

' w t + w x + w xxx + Xw = 

- I Pi 0) - J k(x, y)pi ( y)dy 

< < 

w( 0) = w (L) = w x (L) = 0, 
w t + w x + w xxx + Xw = 0, 
s w x {0) = w(L) = w x (L ) = 0. 

Note that the parameter A is the same for the observer and 
the system itself. 

Given that n, n o are continuous maps, invertible and their 
inverse maps are also continuous, the exponential stability of 
(l36l > would imply the exponential stability of the closed loop 
system and therefore the proof of Theorem Q] will be ended. 

In order to prove the exponential stability of (l36l >. we use a 
Lyapunov argument. Let us consider the following function. 



V(t) = V 1 (t) + V 2 (t) + V 3 (t), 

(37) 

where 

Vi = \w(t,x)\ 2 dx 

(38) 

v 2 (t) = j J \w(t,x)\ 2 dx 

(39) 

V 3 (f) = y J \w t (t,x)\ 2 dx 

(40) 

with A, B to be chosen later. 



di(V 2 (t) + V 3 (t)) < V 2 (t) + V 3 (t) < d 2 (V 2 (t) + V 3 (t)). 


o 


Taking the time derivative of the function V(t), we get 
after some computations that 


V\(t) =A / wt(t,x)w(t,x)dx 
Jo 

r L 


< (-AX + D 2 


\w(t,x)\ 2 dx + A 2 \w xx (t,L )\ 2 


— 2^ — A + 4” Al 2 |rZi ;E3 ;(f, L)\~ 

where D := max x€ [ 0i i,] |pi(x) - k(x,y)p 1 (y')dyj. 

Prom the same computation as in ©, we see that 


V 2 (t) <-2AV 2 (t). 

Moreover, thanks to the regularity H 3 (0,L), the same 
computation can be applied to z = w t (see the proof of 
Lemma Q} to obtain 


V 3 (t) < — 2AV3 (f) 

Thus, we get: 

V(t) ^2^ — A + Vi(t) + A 2 \w xx (t , L)| 2 

-2XV 2 {t)-2XV 3 {t). 

We need to find an upper bound for \wxx(t, L)\ 2 . We 
multiply 

f W t + Wx + w xxx + Aw = 0, 

< (42) 

I w)( 0 ) = w{L) = w x {L) = 0 , 

by xw xx and after some computations we get 


\w xx (t, L)f < \\dj xx \\l 2 (o,l) 

+ ^2A+ —^ ||ujx||l2( 0l i) 

+ ^ll^t II 1,2(0,L) 

and finally the existence of a,b> 0 such that 

\w X x(t, L )\ 2 < a||u}|| 2 2( 0ji ) + 6||ui t |||2( 0ii ) (44) 

Remark 2. Here, we have used that the norm ||/||rf3(o,L) 
and the norm ||/||z, 2 (o,z,) + \\fxxx\\L 2 (o,L) are equivalent. See 
also Remark\J\ o 




We use the latter inequality to write: 

V{t) <2^ - A + Vi(t) + 2a—Vz (f) 

A 2 

+ 2 b—V 3 (t) - 2A V 2 {t) - 2A V 3 (t). 

JD 

Therefore, 

/ n 2 \ / nA 2 \ 

V(t) <2 ( - A + —) Vi (f) + 2 ( — A + —) V 2 (t) 

2(-X+ b -^)Vs(t). 

In this way, by tuning A, B large enough, we get for any 
e > 0 that 

V(t) < 2(-A + e)v(i), 

which gives an exponential stability with decay rate as close 
to A as we want. The rapid stabilization is achieved because 
the parameter A can be chosen as large as desired. 

It concludes the proof of the stability of the closed loop 
system with the output feedback control law depending on 
a boundary measurement of the state. 

VII. Conclusions 

In this paper we have designed an output feedback con¬ 
troller for the linear Korteweg-de Vries equation posed on 
a bounded interval. This boundary feedback control acts on 
the Dirichlet boundary condition on the left endpoint and 
exponentially stabilizes the system by using the boundary 
measurement y(t) = u xx (t,L). Because of this choice of 
y(t), we have to work on a more regular framework, given 
by H 3 (0,L) as the state space. The backstepping method 
is applied together a classical observer in order to build the 
control. 

This work is the first one addressing the output feedback 
problem for the Korteweg-de Vries equation and opens 
several possible extensions. Different location for the control 
or different boundary conditions can be considered as in 
[2], [11]. The choice of the measurement is particularly 
interesting. For instance, it would be nice to deal with 
the collocated case, which could be harder than the non- 
collocated problem considered in this paper. Of course, the 
nonlinear case also appears as a natural next step in the 
study of this control system. In that context, even performing 
simulations of the closed loop system may be an interesting 
challenge. 
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